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1. Introduction 

The reggeized gluon states, also called Reggeons, appear in scattering 
processes of Quantum Chromodynamics (QCD) in the Regge limit where 
the square of the total energy s is large while the transfer of four-momentum 
t is low and fixed. In this limit the leading contribution to the scattering 
amplitude of hadrons is dominated by the exchange of intermediate particles, 
Reggeons, which are the compound states of gluons ^ 121 El E] ■ 

Even in the simplifying Regge limit in the generalized leading logarithm 
approximation [HI E] this problem is technically very complicated due to 
the non-abelian structure of QCD. Many Reggeon wave-functions are the 
eigenstates of a Hamiltonian which is equivalent in the 't Hooft's multi- 
colour limit ^01 ^2 ^] to the Hamiltonian for the non-compact Heisenberg 
SL(2, C) spin magnet. Moreover, the multi-Reggeon system with particles 
is completely solvable ^SIE]- Thus, it possesses a complete set of integrals 
of motion, so called conformal charges {q2,Q2iQ3iQ3T ■ ■ ^QNjQn) that com- 
mute with the Hamiltonian. The eigenvalues of the SL(2,C) Hamiltonian 
are also called the energies of the Reggeons. The Schrodinger equation for 
the lowest non-trivial case, i.e. for = 2 Reggeons, was formulated and 
solved by Balitsky, Fadin, Kuraev and Lipatov ^SlEl^- They calculated 
the energy of the Pomeron state with N = 2 Reggeons. An integral equa- 
tion for three and more Reggeons was formulated in Refs. [3 [HIE] in 1980. 
However, it took almost twenty years to obtain the solution for N = 3, which 
corresponds also to the QCD odderon ^SlElEnj- Finally, the solutions for 
higher A^ = 4, . . . , 8 have been found recently in a series of papers |2n 1221123] 
written in collaboration with S.E. Derkachov, CP. Korchemsky and A.N. 
Manashov. Similar approach was taken by de Vega and Lipatov OK]. 

Description of scattering amplitudes in terms of the reggeized gluon 
states is most frequently used for the elastic scattering amplitude of two 
heavy hadrons whose masses are comparable. Here the amplitude is equal 
to a sum over the Regge poles. In the Regge limit, s — > oo and t = const 
they give behaviour like s"^*^ where a(0) is called the intercept and its value 
is close to 1. On the other hand, the intercept is related to the minimum of 
the Reggeon energy defined by the SL(2,C) Hamiltonian. Thus, evaluating 
the spectrum of the SL(2, C) XXX Heisenberg model we can calculate the 
behaviour of the hadron scattering amplitudes. 

Three reggeized gluon states in i— channel correspond either to Pomeron 
exchanges with parity C = or to odderon states with C = —1. The 
Janik-Wosiek solutions include both of them. Using the duality symmetry as 
well as momentum representation Bartels, Lipatov and Vacca found another 
solution with qs = and a(0) = 1 [SEl Ell EHl EHl EHj • On the other hand 
there was shown in Refs. ^Hl that in order to solve the problem and 
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calculate the intercept one may use Q— Baxter method (3T|. It was applied 
to three Reggeon problem in Refs. |2l.,.22j .23J and Refs. (23123; however only 
in the former ones there was a complete agreement with the Janik-Wosiek 
method. The intercept for the three-Reggeon function was also calculated 
directly using variational method in Ref. ^22J- The results agree with values 
obtained by the exact Q— Baxter method. 

In this work we present the generalized Janik-Wosiek method of con- 
structing the Reggeon eigenstates. Moreover, we calculate the reach spec- 
trum of the the three Reggeon energy and the spectrum of the conformal 
charges {q-i-,q^} for = 0, 1, 2, 3. At the end we evaluate corrections to the 
WKB approximation 521|. 

In Section 2 we discuss the current state of knowledge. Next we describe 
construction of the Reggeon eigenfunctions that consists in solving the the 
differential eigenequations of the conformal charges. We systematize the 
knowledge about the ansatzes for eigenstates of {^3,^3} for three reggeized 
gluon, extend the calculations to an arbitrary complex spin s and derive 
differential eigenequations for the conformal charges §3. Moreover, we show 
solutions to the differential eigenequations with = 3 and s = and resum 
obtained series solutions for the ^3 = case. We construct an exact solution 
to the 53— eigenequation, which may be solved by a series method, and find 
the quantization conditions for {(?3,^3} which come from single- valuedness of 
the Reggeon wave- function. The numerical results of this method for = 3 
are presented in Section 4. They agree with solutions found using Q— Baxter 
method (221122]. We present quantized values of ^3 for different Lorentz spins 
n/t = 0, . . . , 3 as well as corrections to the WKB approximation. At the end 
we make final conclusions. 

2. Hamiltonian method 

2.1. The formalism 
In the high energy Regge limit 

s — > 00 t = const (2.1) 

the contribution to the scattering amplitude of two hadrons A and B from 
an exchange of three reggeized gluons may be written as 

AN=z{s,t) = s J d2zoe^^^o-P(|.A(fo)|e-^^^^^/^(929|a2)"'|$B(0)), (2.2) 

where dk = d/dzk, the rapidity Y = Ins/so and {did2d^^ are gluon 
propagators. Reggeized gluon states, which are also called Reggeons, are 
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effective particles. They interact with each other and propagate in the 
channel. The Hamiltonian TL^ is given [3 El IS as a sum of = 3 
BFKL kernels [inilinilil- The wave-functions \^a{B){zq)) = ^A{B){zi - zq) 
describe the coupling of three Reggeons to the scattered particles. The zq— 
integration fixes the momentum transfer t = —pj^ whereas the operators 
1/9^ stand for two-dimensional transverse propagators. 
Defining the functions ^{z) as 

$(z) = {-ifd,,d,,d,.M^) (2.3) 

the scalar product in the amplitude H2.2|l can be rewritten as 

= («>A(zo)|e-"=^«3^^°'"''/4|cj,^(o)) . (2.4) 

The Hamiltonians, TL^ and TL^, are invariant under the coordinate transfor- 
mation of the SL(2,C) group 

^, ^ azk + h _, ^ azk + b 

czk + d' czk + d ' ' 

with k = 1, 2, 3 while ad — be = ad — be = 1 and they are related to each 
other as _ _ 

^(.=o,s=i) ^ (^Q^Q^Q^-^ ^3 (818283) . (2.6) 

Now, one may associate with each particle generators of transformation 
H2.5p This generators are a pair of mutually commuting holomorphic 

and anti-holomorphic spin operators, Sa and 5„ . They satisfy the stan- 



dard commutation relations 



(k) ^(n) 

/3 



ieafSyS''"' S^'^^ and similarly for 

, 

. The generators act on the quantum space of the k-th particle, V^'^'''^' 

as differential operators 

= Zk8,, + s , 5i') = -8,, , 5f = zl8,, + 2szk , 
So = Zk8,,+s, S^^^ = -8^,, =zl8^,+2jzk, (2.7) 

where S^^ = sl'^^ ± ^52*^'' while the complex parameters, s and s, are called 
the complex spins. Thus, the Casimir operator reads 

))2 = (5i'))2 + (Sf + Si'^sf )/2 = s{s - 1) (2.8) 

j=0 
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and similarly for the anti-holomorphic operator {S )^. The eigenstates of 
the SL(2,C) invariant system transform as ^34| i35i| 

^{zk,zk) ^'{zk,zk) = {czk + d)-^'iczk + d)-2^^'(4,4). (2.9) 

Due to the invariance H2.5p of the system we can rewrite the Hamiltonian 

as 

7^3 = ^3+^3, [Hs,H3] = (2.10) 
in terms of the conformal spins H2.7|l 

TV 3 

= Y,H{Jk,k+i)^ Hs = Y,H(JkMi)^ (2-11) 

k=l k=l 

where 

H{J)='^{l-J)+i;{J)-2i;{l) (2.12) 

with ipi^x) = dlnT{x)/dx being the Euler function and J^^i = J^^i. Here 
operators, Jk,k+i and J^^k+i, are defined through the Casimir operators for 
the sum of the spins of the neighbouring Reggeons 

Jk,k+i{Jk,k+i - 1) = {S^"^ + S(^+i))2 (2.13) 

with Sa^ = Sa \ and Jk,k+i is defined similarly. 

In statistical physics H2.1H1 is called the Hamiltonian of the non-compact 
SL(2, C) XXX Heisenberg homogeneous spin magnets. It describes the 
nearest neighbour interaction between three non-compact SL(2,C) spins at- 
tached to the particles with periodic boundary conditions. 

In QCD values of (s, s) depend on a chosen scalar product in the space of 
the wave-functions H2.9|l and they are usually equal to (0, 1) or (0, 0) |33ll2^ . 

In order to find the high energy behaviour of the scattering amplitude 
we have to solve the Schrodinger equation 

ni'='^''=^^^izi,Z2, i*3) = EN^izi, Z2, fa) (2.14) 

with the eigenstate "^{zi, Z2, z^) being single-valued function on the planes 
= (•Zfcj^A;); normalizable with respect to the SL(2,C) invariant scalar 
product 

II^IP = (^1^) = J d^Zld^Z2d^Z3\^izi,Z2,Z3)\\ (2-15) 

where d'^Zi = dxidyi = dzidzi/2 with Zi = z* . 

Prom the point of view of the SL(2,C) spin chain 

7^3 = 7^(^=0.^=0) ^ ^2.16) 
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Indeed, the transformation Sa — > {018283) Sa {818283) ^ maps the 



SL(2, C) generators of the spin s = into those with the spin s = 1. 

One concludes from H2.15|l that the Hamiltonian -}{^~^'^~^^ jg advanta- 
geous with respect to Tis = as it has the quantum numbers of 
the principal series representation of the SL(2, C) group. 

However, one can also use or even 

the factor {818283)^^ or {818283)^^ has to be included in the scalar product, 
i.e. for (s = l,s = 1) we have 

\m^ = J (fzi(fz2(f Z3\ {818283)'' ^{zi,Z2,Z3)f. (2.17) 

Here the scalar product is no longer in the principal series representation of 
the SL(2,C) group. All these Hamiltonians with the corresponding scalar 
products are equivalent up to the zero modes of the operators {818283)^^ 
and {818283)^ . 

A part of the amplitude (|2.2|1 for = 3 reggeized gluons describes 
the leading contribution of the states with parity C = — 1, odderon, and 
subleading contribution to the C = +1 states related to the Pomeron. Both 
contributions are of the same order. 

Instead solving the Hamiltonian eigenproblem (|2.14p we can solve eigen- 
problems for conformal charges 

N / 2 

*2>n=l \ji=0 
N 

= E {{z^,^,f^'-'^^z^^8z^^{z.,,,r + 2s{s-l)), 

i2>n=i 
N 

/: _ 9 ^ cin) oi-i-i) q(i3) 

13 — ^ ^n«2«3'-'o '-'1 '-'2 

*1,*2,«3 = 1 

N 

= ^ ^ ^ ^^jii2^j2»3^«3«l^^i3^2i2^-^n ~'~ '^^*l*2(^«2«3 ~ ^«3«l)^2;i2'9zij^ 

j3>«2>«l=l 

+ SZijis (^isii — Zi^i^8 z^Ji + SZi^i-^{zi^i-^ — Zi-^i2)8z^,^8zi^ 



2'S Zi^i28z^^ 2s Zi2i-)^8z^_^ 2s Zi^ij^8z^^j , (2.18) 

where Zij = Zi — Zj. Similar relations hold for the anti-holomorphic sector. 
This gives for the SL(2,C) spin s = 

^3 = -iZi2Z23Z3l8z^8z28zs (2.19) 
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and for s = 1 

93 = -idz^dz2dz-iZi2Z22.z:ii . (2.20) 
The eigenvalues of the quadratic conformal charge read: 

^2 = - 1) + 3s(s - 1) q2 = -h(h - 1) + 3s{s - 1) (2.21) 

where conformal weights satisfy 

h = — \-ii^h, h= — \-iUh, (2.22) 

while the eigenvalues of the cubic operator = q^. The parameter has 
the meaning of the two-dimensional Lorentz spin of the particle, whereas 
Uh defines its scaling dimension. They define the transformation law of the 
Hamiltonian eigenstates 

*(ii'o', ^2'0', %o') = {czq + df^iczo + df'' 

X (Yliczk + df^^iczk + df'"] ^-(^10,^20,^30) (2.23) 



with 



^j^{zi ,Z2,Z3) = J d'^Zoe'^'^-P^izi - Zo, Z2-ZQ,Z3- zo) . (2.24) 

The eigenstates ^'(zio, ^20, -?3o) belonging to V^'^''^^ are labelled by the 
centre-of-mass coordinate zq and can be chosen to have the SL(2,C) trans- 
formation properties with zq and zq transforming in the same way as and 

zk, (EISl. 

Conformal charges commute also with cyclic particle permutation oper- 
ator 

r^g,g{zi,Z2, Z3 ) ='^-,,5(^2 , Z3, Zi) = 6^^^ ("'^^ ^-,,5(^1 , £2 , ^s) , (2.25) 

where the conformal charges are denoted by g = {q2, qs) and q = (92,93)- 

The phase Osiq) which is connected with eigenvalues of P is called the 
quasimomentum. It takes the following values [2^ 



e3iq,q) = 27T^, for fc = 0,l,...,2. (2.26) 
The eigenstates of the conformal charges qk diagonalize Ti. and F. 
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Let us also define another operator 

M^'±(fi,f2,?3) = ^^(^3,?2,?l) = ±^^{zi,Z2,Zs), (2.27) 

so called mirror permutation operator which has two eigenvalues, ±1. The 
M operator commutes with the Hamiltonian but it does not commute with 
q2 and qs. 

The cyclic and mirror permutation symmetries come from the Bose sym- 
metry. Physical states should possess both symmetries. 

It turns out that adding colour factor corresponding to the antisymmetric 
constant fabc for the odd-mirror states and the symmetric one date for the 
even-mirror states, we are able to restore Bose symmetry. The tensor fabc 
corresponds to the Pomeron states while dabc is related to the odderon states. 
Therefore, in order to check a C— parity of a given state we need to study its 
parity under the mirror permutation M. The states ^ satisfying M^' = — ^' 
are the Pomeron states whereas states for which M ^ = -t-^* are the odderon 
states. 



3. Various ansatzes 

Using the eigenequation for q2 we get Lipatov's ansatz [2^1 for holomor- 
phic part of the eigenstates: 

1 f ^31 \ '^"^s 
^(^10,2^20,^30) = 7 ( —] F{x), (3.1) 

where x = fj^f|^- If we substitute F{x) = G{x) ( _J j we get an- 
other ansatz 

^'(210,220,2:30) = 7 ^7(7 — Nf) G{x). (3.2) 

{zi2Z23Z3ir V (210)^220)^230)^ / 

As we can see we obtained a totally symmetric ansatz which is equivalent 
to the original one. 

Both ansatzes have advantages and disadvantages. The symmetric one 
H3.2p is appropriate if we want to deal with the particle symmetries. The orig- 
inal one H3.H1 has a simpler structure, it contains powers of h (not h/3), so 
we can use it when we want to construct proper single- valuedness conditions. 
One can easily notice that z'^z^ is single- valued because h — h = rih ^ Z,. 
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3.1. Derivation of differential equations for N = 3 
Prom the 53— eigenequation we get the third order differential equation 

iqsFix) = (3s - h){h - 1 - s)(l - 2x)F{x) + {{{h - 2){h - l){x - l)x 
+ 5^(2 + 11(3; - 1)3;) + s(2 - 2h{l - 2xf + ll(x - l)x))F'{x) 

+ {2 + h- 3s)(l - x)x{2x - l)F"{x) + {x - lfx'^F^''\x) , (3.3) 

where 53 is a complex eigenvalue of 53. Prom the QCD point of view the 
most interesting cases are for s = 0: 

iqsFix) = {h- l){h - 2)x{x - l)F'{x) 

+ {h- 2){x - l)x(l - 2x)F"{x) + x'^{x - ifp'^^Xx) (3.4) 

and for s = 1 

iq^F{x) = {h- 3){h - 2){2x - l)F{x) + 

+ ((4 -2h-{h- 8){h - 3)x + {h- 8){h - 3)x^)F'{x) 

+ (5 - h){x - l)x{2x - l)F"{x) + {x- lfx^F^^\x) . (3.5) 

The first such solution was derived and found numerically in [TH] , 

We solve the above equations by the series method choosing represen- 
tation {s, s) = (0, 0) which gives the same equation in holomorphic and 
antiholomorphic sectors. The differential equation H3.4|l has three regular 
singular points at x = 0, x = 1 and x = 00. To generate solutions around 
other singular points we exchange variables. For the case a; = 1 we can use 
a substitution x = 1 — y: 

iq-iG{v) = {h-l){h-2)y{l-y)G'{y) 

+ {h- 2)(1 - y)y{l - 2y)G"{y) - y\y - \fG^'-\y). (3.6) 

with F(l — y) = G{y) and for x = 00 we have x = 1/y: 
q^G{y)=i{{y-l){h + l){h-2y)G'{y) 

+{y - i)y(2(/i + i)-{h + A)y)G"{y) + y{y - ifG^^'Xy)) . (3.7) 

where F{l/y) = G{y). 

In Ref. PH] Eq. H3.4|l was solved for the case with h = h = 1/2 and 
Refe] = 0. Next in Refs. jSEl EH EHl states with Uh = 0,2 and i^/i G M 
as well as = and ii'h G ^ were found. In this work we extend these 
calculations to = 0, 1,2,3 and & To get a broader perspective on 
the problem the reader is referred to Ref. |39| . 
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3.2. Wave-function for s = s = around x = x = 

In order to obtain the full-complex solution for three- Reggeon state con- 
taining the holomorphic and anti-holomorphic parts we have to glue together 
solutions in these sectors: 



=%({za)^-AW 93,^3) (3.8) 

where we use a (3 x 3) mixing-matrix, A^g^, jIHl which does not depend 
on particle coordinates but only on q = {q2,Q3}- From the QCD point of 
view we have two possibilities of gluing solutions: (s = 0,s = 0) [21] and 
(s = 0, s = 1) [3S]- These two cases are equivalent except for the zero modes 
of the highest conformal charge qn- Let us consider the first case, s = s = 0. 

The antiholomorphic and holomorphic conformal charges are related by 
conditions h = 1 — h* and q^. = q^.. The wave- function has to be single- 
valued. This condition defines the structure of the mixing-matrix. 

For /i Z and / we have solutions of the following type 

oo 

ui{x) = x^y^^an^nx"- , 

n=0 

oo 

n=0 

oo oo 

"Log(x) (3.9) 

n=0 n=0 



and similarly in the anti-holomorphic sector. The recurrence relations for 
the ttn^n are given in Appendix C. One can notice that x°'x^ is single- valued 
only if a — 6 e Z. Moreover we have also terms with Log(rc) which have to 
give in a sum Log(a;x). So in this case we have a mixing matrix of the form 



A^''\h,h,qs,q^) 



a 
/? 7 
7 



(3.10) 



where a, /3, 7 are arbitrary. In the above matrix we have A12 = A13 = 
A21 = A31 = in order to eliminate multi-valuedness coming from the 
power-terms, A23 = A32 to obtain single-valuedness in Log(x)— terms and 
^33 = because the term Log(a;)Log(x) is not single- valued on the x— plane. 
In the case of 93 = and h {0, 1} we don't have any Log(x)— terms so 
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the mixing matrix looks like ^ 

^W(/l,7I,g3 = 0,^3 = 0) 



/3 







a p 

£ 7 



(3.11) 



For / and h £ Z we have a solution with only integer powers of 
X and solutions without a Log, with one-Log and with a double-Log. The 
structure of the matrix is 



a /3 7 
/? 27 
7 



(3.12) 



In the last case for qs = and h £ {0, 1} we have only integer powers of 
X and the third solution with one-Log term. The matrix has a form 



G {0,1}, he {0,1}, qs = 0,^3 = 0) 



070 

13 p £ 
£ 



(3.13) 



3.3. Wave-function for s = s = around x = x = 1 and x = x = 00 

We construct the wave-function around the other singular point exactly 
in the same way, obtaining matrices A^^ \h,h,q^,q^)^ A^^ \h,h,q-i,q^) 
(around 1) and A^°^ \h,h,qz,q'i) (around 00). 

Thus, we have the wave-functions similar to (|3.8p . For h ^Z, and ^3/0 

get 



ui{x) 

U2{x) 

u^{x) 



(l-x)'^^a„,,,(l-x)", 

n=0 
00 

(l-x)^ J]a„,^2(l-x)", 

n=0 

00 00 
(1 _ J2 bn,r,{l -Xr + {1-X)^Y1 - ^)"L0g(l 



(3.14) 



n=0 



n=0 



and similarly in the anti-holomorphic sector. The mixing matrices take the 
following forms 



A^'\h,h,qs,qs) 



a 

P 7 
7 



(3.15) 



^ Greek variables in each A— matrix have different numerical values 
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A«(/l,/i,g3 = 0,^3 = 0) = 



AW(/i€Z,/iGZ,g3,g3) 



13 











a 


P 





e 


7 


a 




7 


P 


27 





7 









(/i G {0, 1}, h G {0, 1}, ^3 = 0, qs = 0) 



070 

(3 p £ 

e 



(3.16) 



(3.17) 



(3.18) 



Similarly, we proceed around x = 00 . For h ^ "L and §3 7^ we have 
solutions of type 



ui{x) = (l/x)0^a„,,,(l/x)", 

n=0 

00 

U2{x) = (l/x)i-'^^an,,,(l/x)", 

n=0 

00 00 

u^{x) = (l/x)-'*^5„,,3X« + (l/x)i-''^a„,,,x«Log(x) 



(3.19) 



n=0 



n=0 



and similarly in the anti-holomorphic region. In this case we have the ma- 
trices 



A^^\h,h,q3 = 0,q^ = 0) 
A^^\heZ,heZ,q:,,qs) 

^(°°) (/i G {0, 1}, ^ G {0, 1}, 53 = 0, ^3 = 0) 



a 















7 







7 









" f3 





" 







Q! 


P 









. 




a 


P 


7 




P 


27 







7 









q; 7 
P p £ 
£ 



(3.20) 



(3.21) 



(3.22) 



(3.23) 
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3.4- Transition matrices between solutions around different poles 

The above solutions around a; = 0, l,oo have a convergence radius equal 
to the difference between the two singular points: the point around which 
the solution is defined and the nearest one of the remaining two. In order to 
define a global solution which is convergent in the entire complex plane we 
have to glue the solutions defined around different singular points. This can 
be done by expanding one set of solutions in terms of the other solutions 
in the overlap region of the two considered solutions. Thus, in the overlap 
region we can define the transition matrices A, F ^ where 



(x,g) = Z^(<7)n«(x,g), 



(3.24) 



Matrices, A and -T, are constructed in terms of the ratios of certain 
determinants 001 ■ For example, to calculate the matrix A we construct a 
Wrohskian 



W 



/(I) 



u {x; q) u 2 {x; 



4^ 



{x;q) 



u'f\x;< 



/(I) 



{x;q) 



4«( 

//(I) 



x;q) 
{x;q) 



(3.25) 



Next we construct determinants Wij which are obtained from W by replacing 
j-th column by the i-th solution around x = 0, i.e. for i = 1 and j = 2 we 
have 
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U 



{x;q3) uf\x;q3) u)^''{x;q3) 



u\'{x;q3) u\'{x;q3) 
u I '{x;q3) u"\ '{x;q3) 

The matrix elements Aij are given by 



u 



(1) 



,//(!) 



{x; qs) 



(3.26) 



A 



W 



(3.27) 



Matrix A does not depend on x, but only on q^. In the similar way we can 
get matrices F and their anti-holomorphic equivalents: A, F . 

Substituting equation l|3.24p into the wave-function (|3.8|1 . one finds the 
following conditions for continuity of the matrix A(g, g): 



A^{q3)A^^\q3,q3)A{q3) = ^(^(gg.ga) , 
T'^{q^)A^^\q.,,q3)F{q3) = A^^) {q^, q^) . 



(3.28) 
(3.29) 



Each Equation, (|3.28|l - (|3.29p . consists of nine equations. Solving them nu- 
merically, we obtain values of parameters a, /3, 7, . . . as well as quantized 



14 



phdnSa printed on 1st February 2008 



values of the conformal charges, qk and qj^. We have verified numerically 
that the spectrum of qk obtained by the above method is equivalent to the 
spectrum obtained using the Baxter Q— operator method which is presented 
in next Section. 

3.5. Additional conditions coming from the particle permutation invariance 

Multi-Reggeon states have additionally the cyclic H2.25p and mirror per- 
mutation (|2.27p symmetries. The conformal charges commute only with P. 
Thus, the eigenstates of q2 and (j2.18p are usually not eigenstates of M, 
so they usually have mixed C-parity. Therefore in order to get solutions of 
a given C-parity we have to impose mirror symmetry. Let us illustrate this 
for the case of 53 = 0. 

For this case we can easily resum the series solutions, see Appendix C. Let 
us take the case for h {0, 1}. The eigenequation for the cyclic permutation 
P with quasimomentum 03((?) gives the following condition 

w'^uf' (/? + j{-xf{-xf + a{x - lY'{x - if 

+pi-xf(x - If + six - lf{-xf) = 
= e'^^^'i^w'^uf' (a + Pi-xf{-xf + -f{x - lf{x - if 

+p{x-lf + eix-lf^ . (3.30) 

Here we have used cyclic transformation laws ()A.6|1 . 

Comparing these two lines we obtain conditions: a = e~*^^*^*)/3, /3 = 
g-i6'3(q)^^ 7 = e~^^^^'^'>a and p = e = 0. One can derive that exp(3i^3(g)) = 1 
so 6*3 ((7) = ^ where k = 0,1,2 (k = for physical states). Thus the 
eigenstate of IP [27j reads 

*(^10,^20,^30) = 

= w''w'' (l + e'^{-xf{-xf + e'^{x-lf{x-lf'^ , (3.31) 

with k = 0,1,2 and where we have omitted the normalization constant. 

Now we can act with a mirror permutation operator on (|2.25p and test 
its eigenequation H2.27|l . Using the mirror transformations ()A.7|1 . similarly 
to (|3.30p . we obtain the following relation 

w^w^i-l)""^ (e'^ + {-xf{-xf + {x - lf{x - if^ = 

= ±w^w^ [l + e'^i-xfi-xf + e'^{x-lf{x-lf'^ . (3.32) 
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Comparing both sides of iHO^ gives (-1)"'^ exp(i^) = ±1, (-1)"'^ = 
±exp(i^) and (-1)"'^ = ±1 where the SL(2,C) Lorentz spins nh = h-h. 
These conditions are consistent with k = 0, |. Only the first case agrees 
with the cychc permutation condition. As we can see for odd we have 
minus sign, so taking into account colour factors ( — 1)^, solution (|;-i.H2|l is 
C— even. For even we have plus sign thus solution is C— odd. The last 
case is unnormalizable because when x — > or x — > 1 it does not vanish so 
the norm, with (s = 0,s = 0), is divergent |27| . 

Using the duality symmetry [2^1 [2H1 IIHI OD]; which corresponds to 
h ^ 1 — h, Bartels, Lipatov and Vacca constructed an eigenstate with ^3 = 
and C = -1 



*(?10,^20,^30) = 

= w^w^x{l - x)x{l - x) (^5(2) (x) - 5(2) (1 -x) + ^5|J<5(^^ (y 



X . 

(3.33) 



This wave-function cannot be constructed using the method described here 
because it contains non-analytical function, the Dirac delta 5^'^\x). 

Now, let us take the second wave- function with five parameters i.e. for 
53 = and h G {0, 1}. For h = 1 the vector of linearly independent solutions 
reads 

u{x) = [1, (-x), (-x)Log(-x) + {x- l)Log(x - 1)] (3.34) 
and for /i = it is 

u{x) = [Log(x - 1), 1, Log(-x)] . (3.35) 
Combining them we obtain 

^-(210, Z20, i'so) = w^w^ (aLog(x -I) +13 + 7(-x)Log(x - I) + p{-x) 

+e{{-x)hog{-x) + (-x)Log(-x) + {x- l)Log(x - 1))) . (3.36) 

Like in the previous case we write the eigenequation for the cyclic permuta- 
tion 

y;'^(=i)^M=o) ((a - 7 - e)Log(x -l) + {p- (3) + a{-x)hog{x - 1) 
+{-(5){-x) + e(-x)Log(-x) + a{x - l)Log(-x) 
-|-e(— x)Log(— x) -|- 7Log(— x) -|- e(x — l)Log(x — 1)) = 

= e'^^'^'^^w (aLog(x - I) + (3 + 7(-x)Log(x - I) + p{-x) 
+e((-x)Log(-x) + (-x)Log(-x) + (x - l)Log(x - 1))) . (3.37) 
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Thus we get conditions: a = e-'^3(9)(ck - 7 - e), /? = e-*^3(<?)(p _ ^ ^ 
6-^^3(9)0,, p = e-»«3(9)(_^)^ = e-'^3('?)(^ _ and e = e'^^^^'ih. We have 
two types of solutions. 

The first one with 6*3(5) = when a = 7 = — e and (3 = p = 0. It has a 
form 

^(^10, Z20, zso) = w {{-x)Log{{-x){-x)) + {x- l)Log((x - l){x - 1))) . 

(3.38) 

We obtained in this way solutions with Log(x)— terms which have not been 
found before. The similar expressions were shown in j21j as asymptotics of 
the g3 eigenfunction. Acting with the mirror permutation operator on IjH.HHp 
we get 

Mw {{-x)Log{{-x){-x)) + {x- l)Log((x - l)(x - 1))) = 

= -w ((-x)Log((-x)(-x)) + {x- l)Log((x - 1)(2; - 1))) . (3.39) 

We obtained minus sign so this state is also symmetric under C— parity {i.e. 
with C = +1). 

Other solutions have 6*3(5) = 27r/3,47r/3. Thus, a = 7 = e = and 
p = — e~*^3('^)/3. The wave-function has a form 

^{zio,Z2o,Z3o)=w{l + xe'''^''y). (3.40) 

These solutions are not eigenstates of the M operator. 

4. Known features of the spectrum 

In this Section we describe the spectra of the conformal charges obtained 
by numerical solutions [^^j. To this end we apply the Baxter operator 
method and the Janik-Wosiek method described in the previous Section. For 
A'^ = 3 reggeized gluons the Hamiltonian method and the Baxter operator 
method give the same results. 

4.1. Advantages of Q— Baxter operator method 

In order to find the spectrum of the conformal charges for higher N and 
the Reggeon energy Ejy one can use the Q— Baxter operator method p?T|l^. 
The Baxter operator Q{u,u) depends on two complex spectral parameters 
u, u and commutes with the integrals of motion, and Qf.. It satisfies the 
Baxter equation: 

iN{u)Q{u,u) = {u + is)^Q{u + i,u) + (n - zs)^Q(n - i,u) , (4.1) 
?Ar(tt)Q(u, u) = {u + is)^Q{u, u + i) + {u- is)^Q{u, u - i) , (4.2) 
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where tj\f{u) is the auxihary transfer matrix 

tNiu) = 2u^ + q2U^-^ + ... + qN 



(4.3) 



and similarly 1^(11). We solve the Baxter equations [22] using the integral 
ansatz which changes Equations 14.11 and 14.21 to differential equations. The 
latter ones can be solved analogically to the Janik-Wosiek method which 
gives Q— Baxter quantization conditions. 

The Reggeon energy is determined by eigenvalues of the Baxter operator, 
Qq^^q^{u,u) and can be written as 



EN{q,q) = -Im^ln 



u'^'^QqAu + ^(1 -s),u + i{l - s)) 



X (5_g__g(u + i{l - s),u + i{l - s)) 



u=0 



with 

±q = (g2,±g3,---,(±)^97v) 
while the quasimomentum Ojsf 



= ^ In 



Qq,qi-is, -is) ' 



(4.4) 



(4.5) 



(4.6) 



4-2. Trajectories 

Solving the Q— Baxter quantization conditions [2^ we obtain continu- 
ous trajectories in the space of conformal charges. They are built of points, 
{q2{i'h)j • • • ) iNi^h)) which depend on a continuous real parameter enter- 
ing g2, I|2.2H1 and (|2.22p . In order to label the trajectories we introduce the 
set of the integers 

£ = {£l,£2,...,^2(7V-2)} (4.7) 

which parameterize one specified point on each trajectory for given h. Spe- 
cific examples in the following sections will further clarify this point. 

Next we calculate the observables along these trajectories, namely the 
energy (|2.14p and the quasimomentum H2.26|l . The quasimomentum is con- 
stant for all points situated on a given trajectory. The minimum of the 
energy, which means the maximal intercept, for almost all trajectories is 
located at Vh = 0. It turns out that the energy behaves around Vh = Q like 

EN{i^h; ^^■"°'^"^) = ii^r""'' + '^Ni^h^ + 0{v^) (4.8) 
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Thus, the ground state along its trajectory is gapless and the leading contri- 
bution to the scattering amplitude around may be rewritten as a series 
in the strong coupling constant: 

A{s,t) ~ -isy2imsf— --rj^U,NmB,N{t), (4.9) 

(as<7Arlns)V^ 

where = OsNc/tt and (,x,N(t) are the impact factors corresponding to 
the overlap between the wave-functions of scattered particle with the wave- 
function of A^— Reggeons, whereas a]y measures the dispersion of the energy 
on the the trajectory around i^h = 0. 

On the other hand, the energy along the trajectories grows with i^h and 
for oo and finally, we have E]\f{uh;£) ~ hivh'. These parts of the 

trajectory give the lowest contribution to the scattering amplitude. 



4-3. Symmetries 

The spectrum of quantized charges • • • , 9Af is degenerate. This degen- 
eration is caused by two symmetries: 

Qk ^ {-ifqk (4.10) 

which comes from invariance of the Hamiltonian under mirror permutations 
of particles, H2.27|l . and 

qk^Qk (4-11) 

which is connected with the symmetry under interchange of the z— and 
^—sectors. Therefore, the four points, i.e. . {qk}, {(— l)'^9fc}) {qk*} and 
{{—l)^qk*} with k = 2, . . . , N , are related and all of them satisfy the quan- 
tization conditions H,S.28|1 . I|,S.29|1 and have the same energy. 



4.4- Descendent states 

Let us first discuss the spectrum along the trajectories with the highest 
conformal charge qn equal zero for arbitrary £ It turns out [23 IlEl 
123 EHj that the wave-functions of these states are built of {N — 1)— particle 
states. Moreover, their energies 0^ are also equal to the energy of the 
ancestor {N — 1)— particle states: 

EN{q2,q3, ...,qN = 0) = EM~i{q2,q?„ ■ ■ -^qN-i) ■ (4.12) 

Thus, we call them the descendent states of the (A^ — 1)— particle states. 

Generally, for odd N, the descendent state with the minimal 

energy Ej\f{q]\f = 0) = has for q2 = 0, {i.e. for h = 0,1), the remaining 
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integrals of motion (73 = . . . = qtv = as well. For h = \ + iuh, i.e. q2 / 0, 
the odd conformal charges q2k+i = with k = 1, . . . , {N — l)/2 while the 
even ones q2k / and depend on Uh. 

On the other hand, for even A^, the eigenstate with the minimal energy 
^7V^~'^^ is the descendent state of the {N — 1)— particle state which has min- 
imal energy with qn-i + 0. Thus, E'jf'^{qj^ = 0) = £:]^!_"i(gAr_i / 0) > 0. 

Studying more thoroughly this problem one can obtain [23 ES] a relation 
between the quasimomentum On of the descendent state and the ancestor 
one Ojsf-i, which takes the following form 

= -e^^^-i = (-1)^+^ . (4.13) 

gjv=o 

Additionally, one can define a linear operator A (JHl EH) that maps the 

subspace vj:^:!^-^ of the (A^—1)— particle ancestor eigenstates with the quasi- 
momentum ^AT-i = vrA'' into the A^— particle descendent states with gjv = 
and On = 7r(A^ + 1) as 

A: _ y^^iv=.(^+i)) . (4.14) 

It turns out that this operator is nilpotent for the eigenstates which form 
trajectories jSH], i.e. A^^ = 0. Thus, the descendent-state trajectory can 
not be ancestor trajectory for (A^+l)— particle states. However, it is possible 
to built a single state jSH] with q2 = Qs = ■ ■ ■ = Qn = 0, i.e. for only one 
point Uh = 0, that has En = and the eigenvalue of Baxter Q— operator 
defined as _ 

Qr°K^)~^> (4-15) 

u 

where a normalization factor was omitted. For = 3 this state corresponds 
to the wave-function defined in H3.38|l . 

Additional examples of the descendent states for = 3 will be described 
later in the next sections. 



5. Numerical results 

5.1. Quantum numbers of the N = 3 states 

In this section we present the spectrum of qs for three reggeized gluons. 
For the first time such solutions for Re [(73] = were obtained in 02]. The 
authors of Ref. 02] used the method of the q^ eigenfunctions described 
in Section 5. Similar quantization condition was also constructed in |26| . 
Solutions with Refgs] 7^ were found in jHHEHl- Moreover, in the latter 
paper the solutions with n/^ / are described. 
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The first solution using the Baxter Q— operator method was described 
in Ref. Later, similar results were also obtained in Ref. j2SI- It turns 

out that the Baxter Q— operator method jSBl for = 3 is equivalent to the 
method of the qs eigenfunctions described in Section 5. For higher A'^ > 3 
only the Baxter Q— operator method was used to find the quantization values 
of qn and the energy. 



5.2. Lattice structure 



Solving the quantization conditions (j3.28|l and l|3.29|l for = 3 and for 
h = ^^2*^ , i.e. with Vh = 0, we reconstruct the full spectrum of q^. It is 

convenient to show the spectrum in terms of q^^^ rather than 53. Since q^^^ 

is a multi- valued function of complex gs, each eigenstate is represented on 

1/3 

the complex —plane by = 3 different points. Thus, the spectrum is 
symmetric under the transformation 



1/3 



1/3 

exp {2Tiik/2,)q^' , where < A; < 3 . 



(5.1) 



Additionally, mirror symmetry (|4.in|l results in a more regular structure. 
For the total SL(2,C) spin of the system h = 1/2, which means Uh = Q 
and Vh = 0, we present the spectrum in Fig. ^ One can easily notice that 



1.5 



0.5 
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. . ^. J j^. ^ 
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-0.5 
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Re[q. 



l/3n 



0.5 1 



1.5 



1/3 

Figure 1. The spectrum of quantized for the system of = 3 particles. On 
the left the total SL(2,C) spin of the system is equal to /i = 5, while on the right 
h — \. Different symbols stand for different quasimomenta 6*3: stars 6*3 = boxes 
6*3 — 47r/3 triangles 63 = 27r/3. 



' ' \'k ' ' ' '.'m' ■ ■ ■ ■ ■ ■ ■ j^.- ■ ■ ■ ■ ■ ■ 
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the spectrum has the structure close to the equilateral triangle lattice of the 
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leading order WKB approximation (|5.8|l . Indeed, apart from a few points 

1/3 

close to the origin, the quantized values of are located almost exactly 
at the vertices of the WKB lattice. The WKB formula (2Q gives 



fe™(^i,^2)]'/' = A;v=3 + i^i2j , (5.2) 

where ii and £2 are integers, such that their sum ii + ^2 is even. Here the 
lattice spacing is denoted by 



A, 





dx 


41/3^ J 


-00 \/l — 



- ^^(^/^) - 0.395175... . (5.3) 



27r 



1/3 

The lattice of extends onto the whole complex plane except the interior 
of the disk with the radius A3: 



iQs'l > As (5.4) 



situated at ^3 = 0. 



1.5 



0.5 



-0.5 



■ * 'A .' , 



'■*/ ':/ *. 



1.5 



0.5 



-0.5 



* . » . .V- 



-0.5 



l/3n 



0.5 



1 



1.5 



Figure 2. The spectrum of quantized q^^^ for the system of TV = 3 particles. On 
the left the total SL(2,C) spin of the system is equal to /i = |, while on the right 
h — 2. Different symbols stand for different quasimomenta 6*3: stars 63 — boxes 
6*3 — 47r/3 triangles 9^ — 27r/3. 

In accordance with (|5.2jl . a pair of integers £1 and ^2 parameterize the 
1/3 . 

quantized values of ■ Going further, one can calculate the quasimomen- 



-0.5 



1/3. 



0.5 



1.5 
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turn as a function of ii and £2- It has a following form 

27r 

^3(4,^2) = Y^i (mod27r). (5.5) 

Thus, as we can see states with the same value of Re[g3 ] have the same 
quasimomentum. In Fig. Q different quasimomenta are distinguished by 
stars, boxes and triangles. 

The same lattice structure is exhibited by the spectra with different n^. 
However, they have different corrections to the leading order WKB approxi- 

1 /3 mm 

mation for ^3 . These spectra are presented in Figs fDI^ The corrections to 
the lattice structure depend on 52 as seen in (|5.8)1 . Since the WKB lattice is 
obtained in the leading order of the expansion for large conformal charges, 

1 \Q2^'^\ ^ ks^^l) the corrections are bigger for lower l^a*^^!. Later, we 
shall discuss some other features of the corrections to the WKB leading order 
approximation. 

As we can see in Figs. I1I2I for /i G Z we have additionally trajectories 
with qs = 0. They are called the descendent states because their spectra are 
related to the spectra for the — 1 = 2 Reggeon states. We discuss this 
point further below. 



5.3. Trajectories in Uh 

In the previous Section we considered the dependence of 53 on rih for 
Vh = 0. However, the spectrum of conformal charges also depends on the 
continuous parameter Vh with h = ^-^j^+ii^h- It turns out that the spectrum 
is built of trajectories parameterized by real parameter I'h- Each trajectory 
crosses one point (star, box, triangle) in Figs. I1I2I An example of three such 
trajectories is presented in Fig. [31 They are numbered by (^1,^2) = (0,2), 
(2, 2) and (4, 2) whereas they quasimomentum ^3(^1, ^2) = 0, Air/S and 27r/2, 
respectively. 

1/3 

The trajectories cumulate at = 0. When we increase u^, tends 
to infinity and the structure of quantized charges starts to be less regular, 

1 /3 m 

especially for trajectories with lower \q^ \. We can see this in Fig. 0]where 
we project trajectories with h = ^ + iu^ on the 1/^ = plane. Here stars 
denote point with = 0, boxes with = 1 and circles Uh = 2. Grey lines 
are drawn to show the projection of the trajectories for intermediate values 

of l^h- 

For n/i 7^ we notice that the spectra start to rotate with i^h- In Fig. 0] 
and we present trajectories with positive = 0, 1,2,3. Due to the sym- 
metry H4.1H1 . which means /i — > 1 — /i* or — > — n/^ with ^3 q^, or 
equivalently, I'h — > —i^h but ^3 — > ^3, the spectrum for the negative Uh is the 
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Figures. The dependence of quantized qsii^h] ^i, ^2) on the total spin h = 1/2 + Wft,. 
Three curves correspond to the trajectories with ^2) = (0, 2) , (2, 2) and (4, 2). 




same as for the positive ones but it rotates in the opposite direction with 

Some of the results presented in this Section were found in earlier works 
[nnilSHI- Trajectories with quasimomentum ^3 = and = were obtained 
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in Refs. jSEllSHl- The case for /i = 2 + ii/^ with quasimomentum 63 = was 
discussed in Ref. [HHI- In this work we additionally analyse the spectra for 
/i = 1 + ih'h and h = 3/2 + ii'h. 

5.4- Energy and dispersion 

For all trajectories in (g2, 93)— space we can calculate the energy of the 
reggeized gluons using Eq. (|4.4p . Example of the energy spectrum for tra- 
jectories from Fig. IHl with h = ^ + ivh is shown in Fig. IHl 

The energy along the trajectories is a continuous gapless function of Vh- 
As we can see the energy grows with rising \i'h\- For n/^ = it has a 
minimum value min^i^E3{i'h; £1, £2) at = 0. In the case 7^ 0, due the 
bending of the trajectories some minima of the energy are moved away from 
z^ft = |3S]- However, the ground state corresponds to the point (s) on the 

plane of (see Fig. ^ closest to the origin. For = 3 the ground state is 
located on the (0, 2)— trajectory at f /i = and n/i = with quasimomentum 
equal ^3 = 0. Due to the symmetry H4.injl it is doubly-degenerated and its 
conformal charge and energy take the following values: 

.^ground ^ ^0.20526 . . . , ^1™"°'^ = 0.98868 .... (5.6) 

In the vicinity of z^/i = the accumulation of the energy levels is described 
by the dispersion parameter (|4.8p (J3 = 0.9082. 
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We show comparison of the WKB result of Eq. (|5.2|1 with the exact 
expressions for 53 at /i = 1/2 in Fig. ^ and Table ^ One can find that 
the expression H5.2|l describes the excited eigenstates with good accuracy. 
In the case where the eigenstates have smaller 53 agreement becomes less 
accurate. Thus, for the ground state with zgs = 0.20526 . . . the accuracy of 
H5.2p is ~ 20%. Obviously, in the region where the WKB expansion is valid, 

i.e. \qy^\ ^ I '72^^15 Eq. (j5.2|l can be systematically improved by including 
subleading WKB corrections. 



(^1,^2) 


/"„ exact \ 1/3 




-£^3/4 


(0,2) 


0.590 i 


0.684 i 


-0.2472 


(2,2) 


0.358 + 0.621 i 


0.395 + 0.684 i 


-0.6910 


(4,2) 


0.749 + 0.649 i 


0.790 + 0.684 i 


-1.7080 


(6,2) 


1.150 + 0.664i 


1.186 + 0.684Z 


-2.5847 


(8,2) 


1.551 +0.672 i 


1.581 + 0.684 i 


-3.3073 


(10,2) 


1.951 +0.676 i 


1.976 + 0.684 i 


-3.9071 



1 /3 

Table 1. Comparison of the exact spectrum of ^3 at ft = 1/2 with the approximate 
WKB expression Ij5.2ll . The last line defines the corresponding energy i?3(0; £1, £2). 
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5.5. Descendent states for N = 3 

One also can notice in Figs. 0] and 13 that for odd Uh we have states 
with g3 = 0. For Uh = and = it has the energy £"3 = 0, so it 
is lower than (|5.6|1 . These states are descendants of the states with two 
Reggeons (2^1 EH [2H1 1211 OHI ■ We constructed them in Section 3 using the ^3 
eigenfunction method. The wave-functions of these states are described by 
HH.H1|1 . (|H.H8|I and HH.4()jl . These states have the same properties and 

the energy as the corresponding states with — 1 = 2 particles, £^3(92 1 Q3 = 
0) =£2(92), with da [El: 



£2(92) = 4 Remi -h)+ ij{h) - 2^(1)] 



8 Re 



, (5.7) 



where -0(3;) = '^^^^i^) and 52 = —h{h — l). Moreover, their wave-functions 
are built of the two-Reggeon states [231231 and the quasimomentum 63 = 0. 
Contrary to the states with ^3 7^ 0, the states with 93 = (|,S.,S,Sp couple to a 
point-like hadronic impact factors 221ES1^5) like the one for the 7* rjc 
transition. 



5.6. Corrections to WKB 

The WKB formula for the lattice structure of the conformal charge gs 
was derived in paper (2^. This formula tells us that for 53 ^ 00 



1/3 r3(2/3) , 



2tt 



|Q(n) 



|Q(n) 



2 00 

fc=3 



|Q(n) 



where 



i7r(2A;-l)/3 



k=l 



and /i, I2, n = {ni, . . . jUn} are integers, while the coefficient 

. 3^/3 , 



(5.8) 
(5.9) 

(5.10) 



where star denotes complex conjugation. 

After numerical calculations we have noticed that better agreement with 
the exact results is obtained for 
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In order to show this, we calculated the values of the conformal charge for 
h = ^-^^ for n/i = 0, 1, . . . , 19. We evaluated numerically with Im[g3] = 0, 

i.e. a'[\ and separately with Ref^a] = 0, i.e. a^^ where the superscript (i) 
and (r) refers to imaginary and real parts, respectively. Then we fitted 

expansion coefficients o^^'*'' for large [gal in the range . . . 5000 with high 
numerical precision. In order to save space in Table [21 we present results 
only for Uh = 0,1,2 and 3. 





coef. 


k = 3 


fe = 4 


fc = 5 


fc = 6 


fc = 7 





(r) 
(i) 


0.509799695633 
3.490200304367 
4.000000000000 


-10.065761318 
20.065761318 
10.000000000 


-76.722084 
104.722084 
28.000000 


-1508.927 
-600.068 
-2108.995 


-44580. 
-40411. 
-84991. 


1 


a, 

(r) , (i) 


-2.585231705744 
6.585231705744 
4.000000000000 


-34.383865187 
44.383865187 
10.000000000 


-209.366287 
237.366287 
28.000000 


-1828.065 
340.913 
-1487.152 


-26404. 
-16641. 
-43046. 


2 


''' (r) ''' 

u 

k 

(r) , (i) 

H +H 


1.336317342408 
2.663682657592 
4.000000000000 


1.2054728742 
8.7945271258 
10.0000000000 


-3.08256385 
31.08256385 
28.00000000 


-31.51219 
110.89046 
79.37827 


-196.87 
377.64 
180.77 


3 


(r) 
"■k 

a) ' + a) ' 

k k 


2.250754858908 
1.749245141092 
4.000000000000 


6.7220206127 
3.2779793873 
10.0000000000 


22.76106772 
5.23893228 
28.00000000 


79.325249 
-0.430149 
78.895100 


268.13 
-75.81 
192.32 



Table 2. The fitted coefRcient to the series formula of ql^"^ H5.8|l with rih = 0, 1, 2 
and 4 

Coefficients for A; = 0, . . . , 2 agree with formula (j5.8p . i.e. oq = 1, 
oi = 1 and a2 = —1, but as we previously mentioned in (|5.1()|l and H5.1H1 . 
the expansion parameter b/ |Q(n)|^ is different. This difference comes from 
the fact that the series formula H5.8p with H5.1()|l is derived in the limit 

1 ^ \l2'''^\ ks''^!- Since in IjS.lip the value of ^2 is much bigger then 2/3 
the value of the parameter b from H5.1ip in the above limit goes to HS.lOjl . 
One may suppose therefore that the factor 2/3 in Eq. H5.1ip as subleading 
was omitted in the derivation presented in Ref. |21j . 

Secondly, we see that the coefficients ak with k > 2 start to depend on 

1/3 

Uh. Thus, to describe the behaviour of ^3 properly, we have to introduce 
a second expansion parameter, for example 52- We can also notice that for 

(i) (t) 

A; > 2 for a few first coefficients + £ ^ and this sum does not depend 
on rih. 

Moreover, we can see that after the numerical fitting we obtain two 

different sets of the expansion coefficients, {a^^^} and {a^^^}, defined in l|5.8p . 

1/3 

for real and imaginary , respectively. Thus, in order to describe full- 
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complex values of in terms of the series (|5.8p we have to use both sets 

1/3 

of coefficients, one for real and one for imaginary part of . Alternatively, 
we can perform expansion with two small parameters, i.e. ^72/12(^)1^ 
1/|Q('^)P- Since the leading terms, i.e. with oq, ai and 02, for real and 
imaginary q^^^ are equal and known analytically, good approximation is 
obtained using Eq. H5.8|l with (|5.1ip and neglecting higher order terms with 

aA:>3- 

6. Summary 

In this work we have considered the scattering processes in the Regge 
limit H2.1|l where the compound reggeized gluon states, i.e. Reggeons, prop- 
agate in the t— channel and interact with each other. We have performed 
calculations in the generalized leading logarithm approximation (GLLA) 
[ZlIHlin]; in which a number of Reggeons in the t— channel is constant. We 
attempted to find a scattering amplitude of hadrons with multi-Reggeon 
exchange. However, a structure of reggeized gluon states as well as their 
properties have turned out to be so reach, complicated and interesting that 
in this work we have focused on description of the Reggeon state properties 
as well as on analysing the spectra of the energy and integrals of motion. 
The case for = 2 reggeized gluons was calculated in Ref . ^3 ^1 HI • Thus 
in this work we focused on A^ = 3 Reggeon states and dependence of their 
spectrum on conformal Lorentz spin nh and scaling dimension l+2i/h defined 
in l(T22|l . 

In order to simplify the problem one applies the multi-colour limit 
which makes the A^— Reggeon system itTTnll SL(2,C) symmetric and 
completely integrable. In this limit the equation for the A"— Reggeon wave- 
function takes a form of Schrodinger equation (|2.14|l for the non-compact 
XXX Heisenberg magnet model of SL(2,C) spins s ^lEHlEl- Its Hamil- 
tonian describes the nearest neighbour interaction of the Reggeons |El 
propagating in the two-dimensional transverse-coordinates space H2.5p . The 
system has a hidden cyclic and mirror permutation symmetry (|2.25p . It also 
possesses the set of the (A^ — 1) of integrals of motion, which are eigenvalues 
of conformal charges [S3], Qk and g^, H2.19|l - H2.2()|l . Therefore, the operators 
of conformal charges commute with each other and with the Hamiltonian and 
they possess a common set of the eigenstates. However, for the A" = 3 case 
the multi-colour limit does not change the Hamiltonian and its spectrum. 
We have the two integrals of motion ^2 and and additionally constant 
two-dimensional total momentum p. 

Eigenvalues of the lowest conformal charge, (/2, may be parameterized 
H2.2H1 by the complex spins s H2.7p and the conformal weight h, where h can 
be expressed by the integer Lorentz spin Uh and the real parameter Vh related 
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to the scaling dimension H2.22|l . Solving the eigenequation for q2 we have de- 
rived ansatzes for A^— Reggeons states with an arbitrary number of Reggeons 
as well as arbitrary complex spins s. Since the = 3 Reggeon ansatz 
separates variables, the 53— eigenequation can be rewritten as a differential 
equation of a Fuchsian type with three singular points (|,S.4|1 [THj . We have 
solved this equation by a series method. Gluing solutions for different sin- 
gular points, (|,S.9jl . (|,S.14p and H3.19|l . and taking care for normalization and 
single-valuedness of the Reggeon wave-function (|,S.8|) we have obtained the 
quantization conditions l|3.28|l - HH.29|l for integrals of motion (92, 92; 93) ^3) 
which we have solved numerically (221 EHI ■ For 93 = ^3 = the series solu- 
tions have a simple form and we were able to resum them. Thus, we have 
obtained analytical expressions for the three-Reggeon wave-functions with 

93 = (lioni - dionii . 

In this work we have calculated the behaviour of the 93— spectrum for the 
conformal Lorentz spins Uh = 0,1,2,3 and the scaling dimension 1 -|- 2ii'h. 
Some results for n/i > were presented before in Ref. (21] for = 1 

and Refs. [HHI for = 3. The quantized values of gg for given and 
fixed Uh exhibit the WKB lattice structure (|5.2|1 . which for N=3 takes a 
form an equilateral-triangle lattice (|5.2jl . Figs. 11151 The non-leading WKB 
corrections move the quantized value of 93 away from the lattice and cause 
that that the quantized values of 93 lie outside a disk located around the 
origin of the lattice, i.e. near 53 = 0. However, for odd Uh there exist states 
with 53 = 0. They are called descendent states because their wave-functions 
are effectively built of = 2 Reggeon states. The non-descendent state 
with the lowest energy belongs to n/i = sector and its energy is positive 

l|5.6|l . This state is double-degenerated and it appears to be the nearest one 

1/3 

to the origin of the lattice. However, the ground state for A^ = 3 is 
the descendent one with = 1 and energy £'3 = (|5.7|1 . Having found 
the exact values of 53 we are able to calculate corrections to the WKB 
approximation H5.1H1 . Table El These corrections differ from the corrections 
obtained earlier in Ref. (2^. The difference seems to be caused by using only 
one expansion parameter r/ for two various conformal charges, g2 and qs. The 
obtained corrections are subleading to the WKB approximation [2^ which is 

an expansion for large values of conformal charges, i.e. 1 <C \q2^^\ <C l^s^^l- 

The above calculations are of interest not only for perturbative QCD 
but also to statistical physics as the SL(2,C) non-compact XXX Heisenberg 
spin magnet model ^1 021 EI- This work opens the way for further stud- 
ies related to the multi-Reggeon states as well as calculations of scattering 
amplitude for concrete processes. 
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Appendix A 

Conformal invariants and other variables 

Let us consider a difference of coordinates {zi — Z2). It changes under 
the SL(2,C) transformation (|2.5p as 

« - = (Stt^ - S^) = + ■ <^-') 

One can see that during the transformation ()A.Hl the factors {czi + d)~^ 
appear in front of the difference. These factors also exist in the transfor- 



mation law of the wave- function (|2.23p . In order to cancel these factors one 
can construct a fraction where the same additional factors appear in the 
denominator and numerator of the constructed fraction variable. 
The fraction variable 

(Zi - Z2){Z3 - Zq) _ 

X = = (2122^32:0) (A. 2) 

[Zi - Z0){Z3 - Z2) 

is invariant under the SL(2,C) transformations H2.5|l . There is only one 
independent invariant for four coordinates [34, 48J. One can see that the 
fractions coming from the SL(2,C) transformation ^A.l\ cancel because 

• the variable is a function of the coordinate differences l|A.ip . 

• the coordinates in the numerator and denominator of the variable x 
are the same. 

Therefore, simplifying the partial fractions we can obtain expression like 
{{azi + b){czj +d) — {azj + b){czi + d)) which with making use of ad — cb = 1 
goes to {zi — Zj). As we can see, we have to build SL(2, C)— invariants from 
differences of the coordinates. 

It is easy to see that performing permutations of coordinates we can 
construct six different dependent invariants 

{Z1Z2Z3Z0) = X, {Z3Z2Z1Z0) = 

{Z2Z3Z1Z0) = 1/(1 -x), (2:1^32:2^0) = l-x, (A. 3) 

(2;32i2;22o) = ix-l)/x, {Z2Z1Z3Z0) = x/{x-l). 



Let us take another product of 2, 



{z[ - z'q){z'2 - z'q) {Zi - Zo){z2 - Zo) 

= {czo + dfw . (A.4) 
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One can see that since in the denominator we have two additional zq variables 
after the transformation we obtain a multiplying factor (czo + d)"^ . Similarly 
for w'^ we obtain 

w'^ = {czo + df'w'' . (A.5) 
Now one can compare this transformation to the transformation of the 



SL(2,C) wave-function H2.2. 

For iV = 3 the variables w = -, — ^ n^ and x = h~^'u?~^°l trans- 

(Z3- Z(,)(Z2 -zO} {zi-Zo)(Z3-Z2} 

form under the cycling permutation (|2.25|1 as 

^ ' (x-l) i-x) ^ ^ ' 
w w{x — 1) — > w{—x) — > w 

while under the mirror permutation: 

W ^ W X ^ w 

For higher we have more invariants. All of them can be constructed 
[2ni from 

r=l r=l k=r+l 



Variables Xr are subject to the two conditions: 



(A.8) 

Hxr = (A.9) 



r=l 

and 



AT N 



5^(-ir n ^'^ = 0- (A.io) 



From ()A.9p we have xi = {—l)^/Y[^=2^r- From l|A.inp we can calculate 



Sr=2(_j0^^rifc=r+l 
Y[k=3 
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Interchangeably, one can derive from ()A.9p = (—1)^/11^1^ and next 

from ifxrnii 

Et?(-irn 



^TTv-i = ^^,1 . (A.12) 



Thus, we see that we have for N particles A^ — 2 independent invariants built 
of the particle coordinates Zi. 

Appendix B 

Solutions for N = 3 and s = around x = 0'^ 

The eigenequation for IjH.Hp is a differential equation of the third or- 
der so around each singular point, x = 0,1, oo, it has three independent 
solutions. Around a; = 0^ we have an indicial equation 

(/i-n-r)(r + n- l)(n + r) = 0, (B.l) 

so its solutions are given by ri = /i, r2 = 1 and = 0. As we can see 
we have two cases when /i Z (one solution with Log(x)) and h G Z (one 
solution with Log(x) and one solution with Log^(x)). As we will see below 
we also have to consider separately solutions with qs = 0. 

Using the same method we can also found solutions around x = 1 and 
oo. However, to save the space we do not present them. 

Appendix B. 1 Solutions for (?3 7^ and h ^Z 
In the first case, i.e. for h ^ Z and qs ^ the solutions look as follows 

00 

n=0 
00 

ra=0 

00 00 
mix) = x''^^bn,r3x'' + x''^J2an,r2x"'Log{x), (B.2) 
ra=0 ra=0 

where ao,r is arbitrary (e.g. equal to 1) 

(ig3 - {h-2r){h-r)r) 
[h — I — r)r{l + r) 
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and m = n -\- r 



_ {h-m + l){h -m + 2)(m - 2) 

0'n,r = JT w T\ 0,n-2,r 

[n — m)[m — l)m 

{iq^ -{1 + h- m)(m - l){h - 2(m - 1))) 
(n, — m)(m — l)m 

whereas 6o,r3 = ■^7^flo,r2 and ^i^^a is arbitrary. One can notice that coeffi- 
cient 6o,r3 is well defined only for ^ 0. Moreover, 

2 + (/i-3)/i-ig3 , , 3/1-8 , 6 + h{h-6) 

^^'^^ 2(2^7^^ ^^'^^ + ^(^"^'^^ + 2(2-/.) (^-^^ 

and 



_ (/i + 1 - - m + 2)(m - 2) 

(/i — m)(m — Ijm 

% - (1 + - m)(m - l)(/i - 2(m - 1)) 

7Z \7 i\ t'ri-l rs 

(/i — mj(m — Ijm 

/i^ + h{7 - 4m) + (m - 2)(3m - 4) 

^ 77^ \7 T\ On-3,r3 

(/i — mj(m — Ijm 
/i2 - 6/i(m - 1) + 6(m - 1)^ 2m - Sm^ + /i(2m - 1) 

[n — m)[m — Ijm {n — mj{m — Ijm 

(B.6) 

where m = + n. 



Appendix B.2 Solution with = 



In this case ao,r defined in l|B.2|) is arbitrary (e.^. equal to 1). For the 
first solution with ri = /i we have ai^n = whereas for the third one ai^r^ 
is arbitrary (let us take 0). It turns out that we do not need Log— solutions. 
The second solution is more complicated. One can derive the exact formula 
for 

-X k-h _ T{l-h + n) 
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and performing summations 

oo 

ui{x) = x^'^ ^^a„^r.^x" = x^ao^n 

n=0 

oo c 



r(l -h + n) 



J(l - h)T(n + 2) 

n=0 n=0 ^ / V I ; 



X 



= -ao,r,i((l-x)'^-l), 

oo 

n3(x) = x''^ ^ an^r-gX" = ao,r3 • (B.8) 
n=0 

Gathering solutions ()B.8p we have 

-u(x) = ^ + 5(-x)'' + C(x - 1)'* (B.9) 

where A, B, C are arbitrary. The above solution was presented by Lipatov 
and Vacca in Refs. pill^ . 

Appendix B.3 Solution with ^3 7^ 0, ^2 = and h = 1 

For /i = 0, i.e. 52 = 0, and / we have a different set of solutions. 
Here we have three solutions of the indicial equation ()B.Hl which are integer 
ri = 1, r2 = 1, r3 = so the solutions are given by 

00 

ui(x) = x^'i^a„,,,x", (B.IO) 

n=0 

00 00 
U2{x) = X''^'^bn,r2x"' + x''^'^an,riX"'l^Og{x) , 

n=0 n=0 

00 00 00 

U3(x) = x''^ ^Cn^rgx" + 2x''2 ^6n,r2a;"Log(x) + x''^ ^an,rix"'Log^(x) , 

71=0 n=0 n=0 

where ao,r is arbitrary (e.^. equal 1) 

r + r2(2r-3) 

Ol,r — 9/-. , ^ : — (B-H) 

r^{l + r)-tq3 



and m = n + r 



(m-3)(m-2)2 ^ (m - 2)(m - l)(2m - 3) - % 

(B.12) 



[m — \ym [m — Ij^m 



r 1 
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whereas ^ 

= 2^-iQ3bo,r2 + ao,ri - Sai^^J (B.13) 

(n — l)^(n — 2) {2m^ — 3m'^ + m — iqs) 

w^yn + 1) n^(n + 1) 

l + 6n(n-l) (n-l)(5-3n) _ (2 + 3n) 

+ n2(n + l) ''"-'''■1 + + n2(n + l) ""-"'''i ~ n(n + l)''"'''^ 

(B.14) 

while co.rs = i^ao,n, Ci^^g is arbitrary and 

C2,r-3 = -:^^93Cl,r3 + ^'0,r2 " 56l,r2 + 3ao,r-i - 4ai,ri > (B.15) 



_ (n-2)2(n-3) -iga - 6 + n(13 + n(2n - 9)) 

Cn,r3 — 7 TTt Cr— 2,r3 "I 7 TTn ^Vi— l,r-3 

n(n — 1)"^ n(n — 1)"^ 



_ 2(n-2)(3n-8) ^ 26 - 36n + 12n^ ^ 

n(n — 1)^ n(n — Ij^ 

6n2 - 8n + 2 2(7 - 3n) 6(3 - 2n) 

n{n — ly n{n — ly n[n — ly 

2(2 - 3n) , 

+ 77ian-i,n- B.16 

n(n — \y 

Appendix B.4 Solution with q-i = (Zs = and h = I 

In this case ao,r is arbitrary {e.g. equal to 1). For the first solution with 
ri = {h = 1) we have ai^^i = and for the third one ai^rs with r3 = is 
arbitrary (let us take 0). The second solution is more complicated. We need 
Log— solutions 

oo oo 

U2{x) = x'"' bn,r2x"' + x''' ^ a„,^,a;"Log(x) . (B.17) 

n=0 n=0 
Using recurrence relations with hi^r2 — 5^0, n (and an^n = for n > 0) 

_ (rz-l)^(n-2) (n-l)n(2n-l) 
^"'^^ - v?{n + l) n2(n + l) ^""'■'^^ ^^"^^^ 
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one can derive an exact formula for hn,r2 = ^1,^2 (n+i)n performing 
summations with arbitrary 6o,r2(= 0) we have X^^q ^"■.''2^" — 2^i,r2(a^ + 
Log(l — x) — xLog(l — x)), so that 



00 



Ui{x) = X^^ ^^On,ria^" = OO.ri = XttQ^^i j 
n=0 

U2(x) = ao,ri (a^ + Log(l — x) — xLog(l — x) + xLog(x)) , 

00 

n3(x) = x''=' ^an^r-gx"- = ao,r3 ■ (B.19) 



n=0 



Gathering solutions ()B.19|1 we have 

u(x) = ^ + S(-x) + C((x - l)Log(x - 1) + (-x)Log(-x)) , (B.20) 
where A, B, C are arbitrary. See also Ref. |24j . 



Appendix B.5 Solution with 7^ 0, ^2 = anrf = 
For /i = 0, i.e. ^2 = 0, and arbitrary 53 / we have three solutions 

00 

ni(x) = x''i ^a„,,.ix" (B.21) 

n=0 

00 CX) 

U2(x) = x^^ ^ bn,r2x"' + ^ a„,r-i a;"'Log(x) 

n=0 n=0 

00 00 00 

U3(x) = x''3 ^ c„,r3a;" + 2x''2 ^ 5„,r-2a;'*Log(x) + x*"^ ^ a,i,r'ia;"Log^(x) 

n=0 n=0 n=0 

where ri = 1, r2 = and = 0. Here ao,r is arbitrary {e.g. 1) 



-ig3 + 2r^ 



and m = n + r 



__(m^2)2 , -ig3 + 2(m-l)3 

whereas 60, r2 = ~i^^o,ri and bi^r2 is arbitrary while 

&2,r2 = ^((2 - %)^l,r2 + 6ao,ri " Soi^r-J , (B.24) 
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{n-2)\ -ig3 + 2(n-l)3 

bn,r2- ^2 ^n-2,r,+ ^l^n - 1) ""'^^'^ 

(n - 2) (3;i - 4) 6(n - 1) 2 - 3n 

an-3,ri H 5 0'n-2,ri H 7Tan-l,n • {D.Zb) 

n^(n — 1) w^in — 1) 

The coefficient Co,r3 = i^(2<io,ri +&i,r2) ^-^d ci^^a is arbitrary and 
1 5 

C2,r3 = ^(2 - ^Q'3)ci,r3 + 861,^2 - 462,r2 + 3ao,ri - 2"l-'-i ' (B.26) 



(n-2)\ , -i gs + 2(n - 1)^ 

n2(n — 1 



Cn— 2,r3 "I" 2/ 1^ l,r3 



_ 2(n-2)(3n-4) ^ 12(n- 1) ^ 2(2 - 3ra) ^ 

n^(n — 1) n(n— 1) 

2(5 - 3n) 12 2(1 - 3n) 

H oT TTan-3,ri H 9an-2,ri H oT TTan-l,ri • (B.27) 

n^(ra — 1) n^(ra — 1) 

Appendix B. 6 Solution with h = and (/2 = 93 = 

In this case ao,r- is arbitrary. For the first solution with ri = 1 we have 
recurrence relations with ai^n = ^o,o,ri 



0'2,ri — 2^1,ri 



1 

n - 1)2 2n 



n + 1 n + 1 

Summing series up we derive an exact formula for 



an-2,ri H —ran-i,ri ■ (B.28) 



00 00 ^ 

ui{x) = x''^ ^ fln.na;" = a; ^ a^^ao.n — = -Log(l - x) . (B.29) 

n=0 n=0 

In the second solution, with r2 = 0, we have arbitrary ai^r2 (let us take 
0) and solution is U2ix) = ao,r2- The third one is with Log— solutions 

00 00 
U3{x) = x'^ ^ri,rsx'' + x'^ ^ a„,,,x"Log(x) . (B.30) 
n=0 n=0 



Using recurrence relations with arbitrary 60. 

^2,r3 = 2 ^^''■3 



r3 



0n,r3 — ^^2 "n-2,r3 H ^^2 f'n-l,r3 (^C.ciiJ 



phdnSa printed on 1st February 2008 



39 



one can derive an exact formula for hn^rz = ^i,r-3n ^'^^ performing summa- 
tions with arbitrary 6o,r3(= 0) we have Yl^=o^n,rzX'^ = ~^o,r3Log(l — x). 
All these solutions look like 

oo 

ui{x) = x'^^ ^an.n^;" = — ai^r3Log(l — x) , (B.32) 

n=0 

oo 

n=0 

oo 

U3(x) = x''=' ^ bn^rga^" + «o,r2Log(x) = -6o,r3Log(l - x) + ao,r2Log(x) . 

n=0 



Gathering solutions ()B.32|1 we obtain 

u{x) = A + SLog(-x) + CLog(x - 1) , (B.33) 
where A, B, C are arbitrary. See also Ref. |24j . 



Appendix B. 7 Solution with 7^ and h = 2 
For h = 2 and / we have three solutions 

00 

= x''i^a„,,.ix", (B.34) 

n=0 

00 00 

U2{x) = x^^ ^ bn,r2x"' + ^ a^^j-i a;"'Log(x) , 

n=0 n=0 
00 00 00 

^^3(3;) = x"^^ C^^T-gX" + 2x''^ 6„^r2^"Log(x) + X^^ an,riX'^^Og^{. 



n=0 n=0 71=0 



with ri = 2, r2 = 1, = 0. Here ao,ri is arbitrary (e.^. 1) 



-iqs + 2r{r-l)ir-2) 

fli.n — r-n TT ao,ri (hi.cibj 

r(r^ — Ij 



and m = n + ri 



(m — 3)(m — 4) —iQs + 2(m — l)(m — 2)(m — 3) 

m[m — 1) [m — 2)[m — Ijm 

(B.36) 

whereas 6o,r2 = ~^^o,ri and 61,^2 is arbitrary while 

^2,r2 = ^{-mbi,r2 + 4ao,ri " Hai.n) > (B.37) 
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(n-2)(n-3) ^ , -iqs + 2{n - 2){n - l)n ^ 

^n,r,- n{n + l) ''"-2.-2+ n(n2-l) ''"-^''-^ 

ll + 3(n-4)n 6(n - 2)n + 4 

7 9 T\ On— 3,ri H I 7 o 7T On— 2,ri 

n(n^ — Ij n(n^ — Ij 

1 - 3n2 
+ — rran-i.n • (B.38j 

Moreover, co,r3 = if^^o.ra and c^^-g = _ l(^££|±^l:i:2l _ and C2,r3 is 

arbitrary while 

C3,r3 = ^(-^g3C2,r3 + 26o,r2 + 8&l,r2 " 2262,r2 + 12ao,ri " 12ai,rJ , (B.39) 



(n-3)(n-4) -igg + 2(n - 3)(n - 2)(n - 1) 

''"•'■^ " n(n - 1) + n(n-l)(n-2) "^""''"^ 

2(3n(6 - n) - 26) 4(11 - 12m + 3m2) 

+ n(n-l)(n-2) "-''^'^ + + n(n-l)(n-2) ^""'•''^ 
2(2-6m + 3m2) 6(n - 3) 

" n(n-l)(n-2) "-^'^'^ " n(n - l)(n - 2) "'"-''''^ 

12 6 , , 

+ —/ TT«n.-3,ri Z ^an-2,n • B.40 

n(n — 1) nyn — 2) 

Appendix B.8 Solution with = and h = 2 

In this case ao,r is arbitrary (e.^. equal to 1). For the first solution with 
ri = 2 we have remaining coefficients an>o,ri = 0. The second solution, with 
r2 = 1, has an arbitrary 01,^2 (we can take it as 0), a2,r2 = and third one 
with arbitrary ai^r-3, 02,r3 (we also set them to 0). All these solutions look 
like 

00 

n=0 

00 



U2{x) = X^'^ ^ an.ra^;" = XaQ^r2 ; 
n=0 

00 

U3(x) = x''^ Qn,r3a^" = ao,r3 • (B.41) 



n=0 



Gathering solutions ()B.4H1 we have 

u{x) = A + B{-xf + C{x-lf , (B.42) 
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where A, B, C are arbitrary. As we can see this solution corresponds to the 
solution with qa = and arbitrary h {0, 1}. For other integer h {0, 1} 
we can observe the same correspondence. 
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